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Using a three-dimensional formalism that includes rela- 
tivistic kinematics, the effects of negative-energy states, ap- 
proximate boosts of the two-body system, and current con- 
servation we calculate the electromagnetic form factors of the 
deuteron up to Q 2 — 6 GeV 2 . This is done both in impulse 
approximation and with a p-n'y meson-exchange current in- 
cluded. The experimentally-measured quantities A, B, and 
T20 are calculated over the kinematic range probed in recent 
Jefferson Laboratory experiments. The meson- exchange cur- 
rent provides significant strength in A at large Q 2 , but has 
little impact on B or T20. 



Theoretically the presence of a D-state in the deuteron 
means that three independent form factors can be con- 
structed: the charge, magnetic, and quadrupole deuteron 
form factors, Fc, Fm, and Fq. These are related to the 
Breit frame matrix elements of the deuteron electromag- 
netic current operator A^ in the three deuteron magnetic 
sub-states | + 1), |0), and | — 1) via the formulae: 
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with r] = — <9 2 /(4Mf), and Q 2 the square of the four- 
momentum transfer to the deuteron. Hence three ex- 
perimental quantities are required to disentangle the full 
electromagnetic structure of this nucleus. Two of these — 
the structure functions A and B — can be obtained from 
the electron-deuteron differential cross-section using the 
usual Rosenbluth separation. These are related to Fc, 
Fq, and Fm, as follows: 
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The third observable of choice is T20, the tensor- 
polarization observable for electron-deuteron scattering. 
To obtain T20 electrons are scattered from a polarized 
deuteron target. T20 is the ratio of a combination of dif- 
ferential cross-sections for electron-deuteron scattering in 
the three deuteron magnetic sub-states to the unpolar- 
ized cross-section. It is related to Fc, Fq and Fm, by: 
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Recent experiments at the Thomas Jefferson National 
Accelerator Facility (JLab) have probed the electromag- 
netic form factors of the deuteron at large space-like mo- 
mentum transfers. T20 has been measured at Q 2 up to al- 
most 2 GeV 2 , B out to about 1.3 GeV 2 , and A to Q 2 = 6 
GeV 2 . At these momentum transfers relativistic kine- 
matics and dynamics would appear to be a necessary in- 
gredient of any theoretical description. Hence there has 
been considerable effort invested in the construction of 
relativistic formalisms for the NN bound state. If the 
constituents of the bound state are understood to be the 
neutron and the proton then this approach is a logical ex- 
tension of the standard nonrelativistic treatment of the 
NN system. Furthermore, regardless of the momentum 
transfer involved, it is crucial that the consequences of 
electromagnetic gauge invariance be incorporated in the 
calculation. Minimally this means that the electromag- 
netic current of the deuteron must be conserved. Indeed, 
the derivation of Eqs. M)-(@) assumed that the deuteronic 
current A^ was conserved. 

Naturally, A, B, and T20 can be calculated using a 
non-relativistic NN interaction which is fit to the NN 
scattering data. This approach has met with a significant 
amount of success (see, e.g. Refs. [|l],[§). ^ n ^ ms n °t e we 
report on an attempt to imitate this approach using a 
relativistic formalism. To this end we construct an NN 
interaction, place it in a relativistic scattering equation, 
and fit the parameters of our interaction to NN data. 
The electromagnetic form factors of the deuteron pre- 
dicted by this NN model are then calculated. The ap- 
proach includes relativistic kinematics, negative-energy 
states, boost effects, and relativistic pieces of the elec- 
tromagnetic current explicitly at all stages of the cal- 
culation. This three-dimensional (3D) technique has 
been developed and applied in Refs. |g[J5). Here our fo- 
cus is on elastic electron-deuteron scattering. First, we 
give a brief review of our formalism in which we dis- 
play expressions for the bound-state equation and the 
current matrix element in the case of an instantaneous 
two-body interaction. Second, we discuss the inclusion 
of meson-exchange-current (MEC) contributions, espe- 
cially the /97T7 MEC, which is known to give significant 



contributions to electron-deuteron scattering. Then, we 
present our results for A, B, and T20. While T20 is repro- 
duced quite well there is significant discrepancy between 
our calculation and the experimental data for A and B. 

A number of alternative 3D relativistic treatments of 
deuteron dynamics exist (see, for instance, Refs. P |lo|). 
Of these, the closest to this work is that of Hummel 
and Tjon |7j], although we eliminate some approxima- 
tions made there. Also, we do not include the uia-f MEC 
in our calculation. Our impulse approximation results 
are very similar to those of Ref . (7| . 

Consider the Bethe-Salpeter equation (BSE) for a 
bound-state vertex function, T: 



r = KG r. 



(8) 



Here K is, in principle, the sum of all two-particle- 
irreducible NN — > NN graphs. The NN propagator 
Go is the product of spin-half Fcynman propagators for 
each nucleon: Go = id\d,2- In studies of this equation for 
the deuteron bound state jll]] the kernel K included a 
set of single-boson exchanges — in analogy to many non- 
relativistic potential models — yielding the "ladder" ap- 
proximation. However, it is well known that in such an 
approximation the Bethe-Salpeter equation does not give 
the correct one-body limit [12j. In other words, if we 
consider unequal-mass particles, and take one of them 
to be very heavy, Eq. (||) does not become the Dirac 
equation for the light particle moving in the static field 
of the heavy one. This limit is only properly treated in 
Eq. (||) if the full set of ladder and crossed-ladder graphs 
is taken for K [Q. In Ref. || we attempted to remedy 
this, and showed that the pieces of the graphs which ap- 
pear in K and are responsible for the one-body limit can 
be resummed so that Eq. (pi) becomes: 



T = C/(G + G c )r, 



(9) 



where the precise form of Gc was derived in [|3||5j] . For ex- 
act correspondence between (||) and (||) we should have: 



K = U + UG C K. 



(10) 



At the level of the one-boson-exchange interaction, where 
K and U have only their lowest-order pieces, we see 
that Eq. (||) defines an improved "ladder" Bethe-Salpeter 
equation, which does have the correct one-body limit: 



r = i^(Go + Gc)r. 



(11) 



This equation is still four-dimensional. One straight- 
forward way to reduce it to three dimensions is to follow 
Salpeter |13| and assume an instantaneous interaction, 
i.e. make the replacement: 
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where q = (qo , q) is the four- momentum of the me- 

(2) 
son. Since AT inst depends only on the three-vector q this 



approximation reduces Eq. ( |11| ) to a three-dimensional 
equation: 



T inst = K^ st (G + G c )T h 



(13) 



Here the three-dimensional propagator (Go + Gc) is ob- 
tained by integrating over the time-component of relative 
four- momentum: 
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Hereafter we always denote integration over zeroth com- 
ponents of relative four-momenta by angled brackets. 

Before examining (Go + Gc), consider the more stan- 
dard equal-time Green's function |l3Jl5|]: 
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where A are related to projection operators onto posi- 
tive and negative-energy states of the Dirac equation, E 
is the total energy, and e^ = (pf + m 2 ) 1 / 2 . The propaga- 
tor (Go) is not invertible |lq] , since it has no components 
in the H — and — h sectors. This is related to the lack 
of a correct one-body limit in the ladder BSE. If we had 
applied the 3D reduction ( |l2| ) to Eq. (@) we would have 
obtained the Salpeter equation, with the non-invertible 
(Go) in the intermediate state. However, adding (Gc), 
which comes from resumming pieces of the crossed-ladder 
graphs before reducing to three dimensions, gives a 3D 
AW propagator: 
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with k 2 a parameter that enters through the construction 
of Gc- This three-dimensional propagator was derived 
by Mandelzweig and Wallace with k° equal to the on- 
shell energy of particle two JL7J . With k^ chosen in this 
way (Go + Gc) has the correct one-body limits as either 
particle's mass tends to infinity and also is invertible. 

In this work we consider the interaction of two parti- 
cles of equal mass, and thus choose a different k 2 , namely 
k 2 = (E — ei + t2)/2. This form avoids the appearance 
of unphysical singularities when electron-deuteron scat- 
tering is calculated [||. It yields a two-body propagator: 
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which is consistent with that demanded by low-energy 
theorems for composite spin-half particles in scalar and 
vector fields ]Tq] . Alternatively, comparing the ++ — > 
++ piece of 



K^ st (G + Gc)K,_ 



(2) 



(18) 



with the amplitude obtained at fourth order in the full 
4D field theory we see that the contribution of negative- 
energy states agrees at leading order in 1/M [||. In other 
words, effects such as Fig. [l]are included in a bound-state 
calculation that employs Eq. (OJ). This is true even if 

— (2) 

only the instantaneous ladder kernel K^J t is used, be- 
cause of our careful treatment of the one-body limit. 



FIG. 1. One example of a Z-graph which is included in our 
3D equation (tL3[). 

Of course, the instant approximation which led to 
Eq. (OJ) is uncontrolled. However, one of the central 
results of Ref. H was that similar three-dimensional 
reductions could be implemented in a systematically- 
improvable way. The numerical work of Ref. || then 
showed that the replacement of Eq. (jlj) by the system- 
atic 3D reductions developed in Ref. 0] had little impact 
on the deuteron electromagnetic form factors. Hence in 
this work we will report only on results obtained using 
Eq. (|13|) , even though "better" treatments of the reduc- 
tion to three dimensions are certainly possible. 

It remains to construct a conserved electromagnetic 
deuteron current. We begin by defining A]i , the impulse 
approximation current obtained when the ladder Bcthc- 
Salpeter equation (g|) is solved: 



4 0) 



iT(P 



-Q)d 1 ( Pl )d 2 (p 2 + Q)jlf>d 2 (p 2 ) 
(1^2), 



T(P) 
(19) 



where T is the solution of Eq. (g), P = p\ +p 2 and P + Q 
are the four-momenta of the initial and final states, and 
j/J 1 is the usual one-nucleon current for particle n: 
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(q n is the charge) . Using the Ward-Takahashi identity as- 
sociated with formally-modified but practically-identical 

form of this current "" '•' ' '■• ' <''"' ''''' 

conserved, i.e. 



|19| it is easy to show that Ajj, is 
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However, in this work we did not begin with the ladder 
BSE. Instead we began with the 4D equation (^). Con- 
structing a conserved impulse approximation current for 
the vertex function which is the solution of Eq. (H) is a 
little more involved. In Ref. ^ we showed how to add 
a piece to the current (19) which results in A^ being 
conserved if T is the solution of Eq. (£]): 



A^ = A^+t(P + Q)Gc^(P). 



(22) 



The explicit expression for Gc^ can be found in Ref. ||. 
With this four-dimensional current in hand we may 
make a reduction of it to three dimensions in an anal- 
ogous fashion to the reduction employed for the bound- 
state equation itself. Once again, this reduction can be 
performed in a systematic fashion, but here we need only 
the results for an instantaneous interaction. In that case 
the "instant" current: 



Amst,fj, — 1 inst t^inst^ *■ '"' 



(23) 



is conserved, provided that rj, lst is the solution of 
Eq. (fl3) . The explicit form of the current employed is [p[ : 

G^Jpu P 2 ; P, Q) = % (di(pi) d 2 ( P2 + Q)j^d 2 ( P2 )) 

+ t(d 1 ( Pl ) 4( P2 + Q)j^dUP2)) + (1 - 2). (24) 

(n) 

Here d n is the Dirac propagator for particle n, and jji 
is the one-body current ( |20| ) . Only the 7 M piece of jj 1 in 
relevant for charge conservation, since the piece propor- 
tional to a^ u is automatically conserved. Meanwhile, d c n 
is a one-body Dirac propagator used in Gc(P) to con- 
struct the approximation to the crossed-ladder graphs. 
Correspondingly, d c n appears in Gc{P + Q)- It does not 
equal d^, even if particle n is not the nucleon struck by 
the photon. Finally, 



#4 = W ; -# ; ); 
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with p 2 = (e(p2), P2)- The current defined by Eqs. (gj 
( p5| ) includes not only the effects of time-ordered graphs 
like that on the left-hand side of Fig. g, but also the 
effects of photons coupling to the Z-graph in Fig. pi 




FIG. 2. On the left we depict the positive-energy-state im- 
pulse approximation mechanism for electron-deuteron scat- 
tering. On the right the pn"/ MEC diagram is shown. 

This defines our impulse-approximation current. De- 
tailed results for this current employed with the solu- 
tions of Eq. (OJ) were presented in Ref. || . In this work 
we move beyond the impulse approximation by includ- 
ing the effects of the pnj meson-exchange current. Be- 
cause of the quantum numbers of the deuteron this is 
generally thought to be the lowest-mass mesonic excita- 
tion which makes a contribution to the electromagnetic 
deuteron current. However, note that there are addi- 
tional pionic currents which should be included because 



we have chosen pseudovector pion coupling. In nonrel- 
ativistic approaches pseudovector pion coupling implies 
that there is a "relativistic" two-body contribution to 
the deuteron's charge operator involving a jtt contact 
term M|. In A this effect plays a more significant role 
than the pnj MEC fl21J] . This dynamics appears in the 
current operator of the formalism derived in Ref . M , but 
we have not included it in this calculation. 
The Lagrangian governing the pnj vertex is: 
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This yields the two-body current depicted on the right- 
hand side of Fig. 0. In calculating this diagram the same 
form factors are employed at the wNN and pNN vertices 
as were used in calculating the NN potential. The pw^f 
current is automatically conserved. 

With all the theoretical pieces of the puzzle assembled 
we now calculate electron-deuteron scattering. The ver- 
tex functions employed are the ones calculated with all 
positive and negative-energy states included, as described 
in Ref. ||. If the negative-energy states are dropped 
the interaction is exactly the Bonn-B potential for the 
Thompson equation, as derived and fitted to NN phase 
shifts in Ref. p3| . This model gives a reasonable fit to the 
TV TV data, and good deuteron static properties, although 
it is not as good a fit as some more recent NN poten- 
tials |2|,^3|;g4| . When negative-energy states are included 
the deuteron pole position changes slightly. To compen- 
sate for this we adjust the cr-meson coupling from the 

value of the fit in Ref. ||, f^ = 8.08, to f£ = 8.55. 

The single-nucleon form factor parametrization chosen 
is that of Mergell et al. (£5|. This parametrization is 
based on vector-meson dominance with constraints im- 
posed on the asymptotic shape of F\ and F 2 using argu- 
ments from perturbative QCD. 

The only remaining freedom in the calculation is the 
choice of the parameters governing the pir^f MEC. The 
coupling gpn-y = 0.563 is extracted from the decay p — > 
7T7. The contribution of the MEC to electron-deuteron 
scattering depends crucially on the behavior of the cur- 
rent operator as a function of Q 2 . In the work of Hum- 
mel and Tjon vector-meson dominance was used to ob- 
tain a pnj form factor given solely by the 10 meson: 
FpwyiQ 2 ) = qtz — -■ This same pir^f form factor is also 
employed in the non-relativistic calculations of Refs. [jl],|| . 
Other calculations have used form factors based on quark 
models |8). Such form factors tend to reduce the contri- 
bution of this MEC, which is also very sensitive to the 
cutoff masses in the irNN and pNN vertices. 

In Fig. we present our results for A. The left panel 
shows the results up to Q 2 ~ 2 GeV 2 and the right 
panel gives results and data over the full range of ex- 
perimental Q 2 . The two JLab experiments which have 
produced data for A [A198,Ab98] are denoted by triangles 
and squares respectively. Note that these experiments 
confirm the trend of the SLAC data of Arnold et al. 
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FIG. 3. The deuteron structure function A. The left panel 
is an enlarged version of the right one. In both the dashed 
line is the impulse approximation result, and the solid line 
is the result with the pivy MEC included. The experimental 
data of Refs. [28,29] are denoted by the open circles, while 
that of Refs. [26,27] are represented by triangles and squares. 

We see clearly in the left panel of Fig. that the im- 
pulse approximation underpredicts the A data for Q = 
1 — 2 GeV 2 . Much of this lack of strength is made up for 
by the pir^ MEC, which gives a curve that reproduces the 
data in the region Q 2 = 2 - 4 GeV 2 . However, this MEC 
calculation then oue?predicts the JLab data at large Q 2 . 
Our pirj MEC contribution to A scales as Q~ 20 when Q 
is large, which is consistent with perturbative QCD. The 
failure of our calculation to describe the data for Q 2 > 4 
GeV 2 implies that the pirj mechanism is too strong at 
large Q in the present model. 
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FIG. 4. The deuteron structure function B. The experi- 
mental data are given by the open circles. The dashed and 
solid lines represent impulse approximation and pvvy MEC 
calculations. 

As seen in Fig. the extant experimental data for B 
(which, as yet, includes no JLab data) [£9],[30| are even 
less well-described. Already at Q 2 < 1 GeV 2 there is 
significant disagreement between our calculation and the 
data. This disagreement is worsened by the inclusion of 
the relativistic version of the pivy MEC: an effect also 
seen in the calculations of Ref. M . There is a significant 
variation in results for this observable at Q 2 = 1 — 2 GeV 2 



in non-relativistic formulations of this problem, partly 
because in this regime B is sensitive to currents associ- 
ated with the short-range piece of the NN interaction. 
These currents, as well as other mechanisms not in our 
model, could play a key role in reproducing the B data. 




Q'lfnf*) 

FIG. 5. The tensor-polarization observable T20. The older 
experimental data are denoted by the open circles, the recent 
NIKHEF data of Bouwhuis et al. by diamonds, and the JLab 
data of Beise et al. by squares. The dashed and solid line 
represent impulse approximation and pnj MEC calculations. 

On the other hand, Fig. H shows that the T20 data 
of Refs. |3^,M is well-described by our approach out to 
Q 2 w 2 GeV . This observable is fairly insensitive to 
some of the dynamics which plays a role in A and B (e.g. 
the single-nucleon form factors). However, it is quite sen- 
sitive to relativistic effects (see, for instance flnOl), so it 
is gratifying that our approach reproduces the data, es- 
pecially that of Ref. |B3] at large Q 2 , so well. 
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